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\ ABSTRACT 


\ 
\ 

Road profiles available for digital terrain simulation at 
TARADCOM are cataloged. The terrain simulation program for the 
PACER 100/EAI 580 Hybrid computer system is described. The terrain 
simulation program was developed to better utilize the hybrid 
computing system when it is required to run several vehicles over 
known terrains at different speeds for varying lengths of time. 
The operator can interactively setup the test environment or feed 
in a prepared paper tape. Terrains are stored on disks and vary 
from profiles of actual test courses to synthesized forcing 
functions. Development and operation of the program are described 
and possible future refinements are explored., 
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PREFACE 


This report describes digital terrain simulation at TARADCOM. 
The aim of the simulation is to subject a vehicle to a vibration 
environment similar to what it might experience in the field. This 
is accomplished by placing the vehicle on hydraulic actuators which 
move up or down depending on the electrical voltages applied. Elec- 
trical signals are usually supplied by an analog computer. Digital 
terrain simulation is accomplished using a digita] computer to 


produce voltages in digital-to-analog (D-to-A) components on the 
analog computer. 


Reference 1 describes how several actual terrain profiles were 
digitized for use as input to vehicle simulation studies. These 
profiles and others (about 40 in all) were transferred to disk 
storage for rapid on-jine use by the digital computer. 


The report describes how these terrain profiles were converted 


to wheel displacements and how the wheel displacements are applied 
to a vehicle. 


1. CONVERSION OF TERRAINS 


The terrain elevations are in feet and were measured at one, three, 
six, and twelve-inch intervals. Most terrains had 600 or less data points 
and this was chosen as a practical limit. All elevations are greater than 
zero and some reach a height of almost eight feet. This presented a 
problem because the actuators currently available at TARADCOM only have 
a 12-inch span. It was solved by taking into account the fact that long 
uphill] or downhil] trends provide very little input to the vibration 
response of a vehicle. 


| 


The detrending routine is fairly straightforward and is presented in 
Appendix A. The program takes out long uphill and downhill] trends by 
° drawing a line through the slope which has no data point more than six 
inches away. 


eee . ” 
Figure } 
TERRAIN ABOVE ZERO 
r The slope is then moved down to level, with the starting and 


ending points equal to zero. 


Figure 2 
DETRENDED TERRAIN 


In this way the relative slopes between data points are preserved 
while keeping within the physical limits of the actuators. The only 
discontinuities exist at the pivot points of the detrended segments. 
The subroutine finds the longest segments in order to minimize the 
number of discontinuities. 


Now that the terrains fit an actuator it is necessary to calculate 
the motion of a wheel traversing the data points. The wheel trajectory 
routine described in reference 2 obtains the position of the bottom of a 
rigid wheel traversing a non-deformable terrain. The resulting wheel 
trajectories for each terrain were then dumped to paper tape and loaded 
into separate disk files on a digital/analog hybrid computer. 


2. SIMULATION PROGRAM 


nese act SMA iMinENSe tale ane nnaiasatalremanmeeceee 


The program has basically three phases: 
1. Terrain Input - Subroutine "INPUT" 
2. Bay Setup ~ Subroutine "BAYS" 

3. Terrain Output ~ Subroutine “RUN" 


In the terrain input phase, the operator is asked the names of the 
files on disks containing terrain data. The first record ofthe disk 
file has the number of samples and sample interval in inches followed 
by records containing the scaled fraction data. In scaled fraction data 
there are 15 bits of information plus a sign bit. For example: 


0000000000000000 = 0 volts = 0 feet 
1000000000000000 #~9.9999 volts = -.5 feet 
0177777777221727 = 9.9999 volts = .5 feet 


After all terrains have been input, the operator is asked how many 
bays will be run. The operator must supply the number of axles and their 
dimensions, which terrains wil] be traversed, the angle of traversal, the 
desired speeds, and distances for each bay. 


The subroutine BAYS translates these parameters into control 
variables which are later used by subroutine RUN. To get an idea how 
this is done, picture a jeep mounted on four actuators. If the wheel 
signal were applied to all four actuators simultaneously the jeep would 
have only vertical accelerations being input into it. Now if the back 


he 


wheel signals were delayed a small amount of time there would be angular 
motion between the front and rear of the vehicle (called pitch). When 
wheels on both sides of a vehicle see duplicate paths the vehicle is said 
to be traversing a “washboard"-type terrain. Traversing the washboard 

at an angle causes a time delay from when a wheel on one side of the 
vehicle hits a point to when a wheel on the other side hits it. This 
produces a side-to-side motion called roll. 


BAYS will translate the speed and angle of attack into relative 
delays between the wheels of a vehicle. In this way vertical accelera- 
tion, pitch, and roll can be produced in a vehicle in the laboratory. 


In the output routine, RUN (see figure 1), the data channels are 
all more-or-less independent. The subroutine scans the channels and 
subtracts the time since the last scan from each channel timer. If the 
timer reaches zero, a new data point is output. A channel may change 
to a new terrain only after a complete run of the current terrain. This 
is to insure that the output of one terrain is zero before a new terrain 
is introduced. Channels common to a bay are sequential and the right 
rear wheel is identified as the LAST to finish the phase. This insures 
proper interwheel timing for the next phase. There is a two-second 
pause between phases to allow easy differentiation of terrains. 


It should be noted that this setup allows two vehicles to use the 
same terrain. If one vehicle changes terrains, it snould in no way 
affect another vehicle which might happen to be in the middle of a run. 


Once the digital program is running, the hydraulics should then be 
turned up to whatever level is desired. To shut down, the hydraulic 
controllers should be brought to zero and the program terminated. 
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3. FUTURE PROJECTS 


a. The system currently has eight channels of analog wheel displacement 
data. This allows, for example, a three-axle truck and a one-axle water 
trailer or, say, a four-axle tracked vehicle by itself. If eight more 
channels were added, there would be greater flexibility in simulation 
scheduling. Three vehicles in the three bays would be virtually independent 
of each other and would not have to double up on data channels. 


b. Due to the digital nature of the terrain simulation, the output is 
a Staircase function as shown below in Figure 4. 


Figure 4 
STAIRCASE FUNCTION 
This introduces high-frequency noise and is not very realistic. A 


straightline approximation as in Figure 5 is desired. There are two 
basic approaches to this problem. 


Figure 5 
LINEAR INTERPOLATION 
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1. Construct tinear interpolators as shown below. 


Figure 6 


LINEAR INTERPOLATOR 


The resources exist to construct only four of these and would necessitate 
changing the potentiometer settings for a bay each time it changed terrains. 
The effect of this would be to require all of the bays to run the same 
terrain at the same speed and for the same distance, which would limit the 
flexibility of the system. 


2. As the number of data points is increased (by decreasing the sample 
: interval) the staircase more closely approximates a straight line. All 
ia is penis is to use Capacitance Limiting to smooth out the smal] 
ncrements. 


This approach is currently being pursued but it is limited by the 
digital computer haying only 16K of memory. The addition of 16K would 
allow a greater number of terrains in a test setup. 


E c. Ina system with a large library of terrains such as this, it soon 
becomes desirable to modify or combine terrains for special applications. 
A system for updating or modifying terrains would be a future enhancement. 


4. TERRAINS 
On the following pages you will find plots of terrain profiles with 


i detrended wheel trajectories below them. The actual wheel displacements 
f are tabulated after each plot. 


The RMS value is calculated by the following formula: 


[oN 
RMS #4) 4 yi - ¥) 2 
Vit my ” 


where: 


N = Number of terrain samples 
Y = Average wheel displacement 


Yi = Wheel displacement for sample i 


The date on the plot is when the terrain was converted. The date 
on the tabulation js when the terrain was first stored on disk. 


TABLE 1 


TERRAINS WITH CORRESPONDING FILE NAMES 


TERRAIN FILE NAME 

Ft. Knox Mild RDO1 

Ft. Knox Rocky RDO2 

Synthetic Field RDO3 

Ft. Knox Medium RDO4 

synthetic Rock and Log RDO5 

APG Perryman III RDO6 

APG Belgian Block RDO7 

Six-Inch Sine Wave Course RDOS 

Three-Inch Spaced Bump RDOS 

RMS Equal One RD1O 

Modified RMS RDI2 

: ‘ APG Terrain 10 RDI3 
APG Terrain 9 RD14 

APG Terrain 1) RDI5 

\ APG Terrain 12 RDI 
Houghton Data 3 RDI? 
Houghton Data 4 RD18 
Houghton Data 5 R019 
Houghton Data 6 RD21 
Houghton Data 7 RD22 
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TABLE 1 


TERRAINS WITH CORRESPONDING FILE NAMES (CONT'D) 


TERRAIN 


TARADCOM Wood Course 
Ft. Knox CTA] 

Ft. Knox CTA2 

Ft. Knox STV1 

Ft. Knox STV9 

Ft. Knox STV4 

APG 30 

APG 29 

APG 34 

APG 35 

APG 37 

APG 32 

APG 36 

APG 34 

APG 31 

APG 29 

Profile 4 APG 
Profile 4 Chrysler 
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COMPUTER PROGRAMS 


1. DETR - Detrends terrain and performs wheel trajectory calculations 


(see section 1 and Appendix C). 


This program reads in a terrain from the disk, graphs it, takes out 
long uphill and downhill trends (subroutine DETRZ), digitally rolls a 
wheel over it (subroutine TRAJ) and then graphs the resulting trajectory 
of the bottom of the wheel. The wheel trajectory can then be listed on 
the line printer and punches out to paper tape. 


DETR is only a control program which handles input, calls the sub- 
routines and provides output. 


Il. DETRZ - Detrends a terrain so that it starts and ends at zero and 


Stays within a specified limit. 
The subroutine is supplied four parameters. 
Y - Real array of terrain data 


N - Number of data points 


Go 
m 
La 
4 
> 
‘ 


Distance between data points 


a 
_ 
nn 
~o 
1 


Maximum distance a data point can be above or below zero 
Using the equation for a line 
Yi =mxXitb 


and tue distance of a point (Xi, Yj) from a line 


D = {Yj - (m Xj + b)} 


re? 
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the subroutine uses Yi as one of the endpoints of a line and scans the 
terrain for another endpoint so that D will be less than DISP. Once 

this line is found, the new terrain elevations are calculated as being 

the distances from this line. The subroutine then uses the second end- 
point as a starting endpoint and scans the rest of the terrain for another 
line. The entire terrain is detrended in this way and the subroutine 
exits. 


III. TRAd - Translates terrain displacements into wheel trajectories. 


This subroutine is taken directly from reference 2. If a rigid 
tire is rolled over a non-deformable terrain, the bottom of the tire may 
not meet al) of the data points. 


The subroutine calculates the height of the bottom of the tire over 
each of the terrain data points. A tire of 42" diameter was used for 
all terrains to obtain maximum smoothing. For a discussion of the 
mathematics involved, see reference 2. 


IV. LOAD - Loads a paper tape containing wheel displacemeits to a disk 
(see Appendix C). 


The program expects the tape to have the five character name of the 
file as the first record, followed by the number of data points and sample 
interval in inches ina 215 format as the second record. The data points 
then follow, ten to a record in a 1055 format. The program creates a 
disk file for each terrain. 


VY. SAMP - Sets up and runs a digital terrain simulation 
(see section 2 and Appendix C). 
SAMP is structured as a contro] program which calls the three subroutines 


INPUT, BAYS and RUN. Communication of the control variables is accomplished 
through the two FORTRAN Common's; DATA and PARAM. 


COMMON DATA 


BEGIN (8) Contains the starting locations of up to eight terrains in 
the array SAMPLE. 


SAMPLE (7500) Contains up to 7500 terrain data points. 
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COMMON PARAM 


: TERRA (8,8) TERRA] jy is the terrain number for CHANNEL; during PHASE). 
\ SPACE (8,8) The time between data points in milliseconds. 
NRUN (8,8) The number of times the terrain will be traversed. 


: DELAY (8,8) The number of milliseconds to wait before outputting the 
: first data point. 


NPHASE (8) NPHASE; is the number of terrains Channel I will cyclicly 
execute. 


LAST (8) LAST; * .TRUE. if I is the lest channel for a bay. 


: INIT (8) If LAST] is true then INIT; will be the first channel for 
: that bay. 
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Agar 
Ags 
289 
aus 
aca 
ugar 
agad 
ags7 
ae 34 
aed t 
Red 
Had 
uae 
ages 
aged 
veas 
Oya 
@u47 
QGa45 
Quay 
ages 
Ov? 
AQS2 
@e53 
Oe54 
teks) 
agha 
G27 
ay 5~ 
nr SG 
aeA. 


6 INL 


[nt -s7Gwa4s phADS A TERRAIN FROM TYRE DIBC, 

VHE FkieST 15 bOKDS AkKE & DESCRIPTION OF TRE TERRAIN, 
wOR? 16 18 [HE MUMBER OF POINTS IN THE TERRAIN, THE. 
REMAINING PORTION CP THE FILE CONTAINS FLOACBAE AOENT 
VARTABLES, 


OF Tre TEPRETM, 


orernanancane 


OhOGRam DOFTR e 
INTEGER TITLE (249) 
Ahat x¥(669),¥ (02) 
FIKESSTION ALEDG VANE (IZA) phAF(S) 
EGU LVALENCE (ay NA), CYeACT AD), CTITLED NA), (Np MAC2IDD 
CALL LPLOT£2) 
2 CONTINUE 
eeTVE (1,129 
lv POR PST CMESTER TERRAIN FILE NAMED") 
Rel (1,20) RAF 
gu FORMaY(2A2,A1) 
CaLt ERABE y 
Cae EXEC(14,.10R25,NA, ueeana hee? 
Oki i4waty?) 7 12, 
IF (N ,6T, 670) NE6uM 
Co Ge iuy,h . 
M(J) es @ DELTA ; 
99 CuntiIrug 
TVEma( 20h) 4260128 
CALL CXFO(14,1%28,NACL29), TIER, NaF I) 
Cau LPLOTII,4&, 030732238, 608,,8,) 
Call S¥MA CU, petererd Se NAF» 95,33 
Cate LPLOT (GX, VoNop tee yo? 
CALL LETRZ(Y,Np DELTA, 5) 
COLL TRATION TFINCAC12)95% 9420) 
CALL LPLOV COs Me VaN ele Mar Bed 
Resi. (1,¢@F) NOYES 
IF yupyeS ,~hk,e 2HYE) GN TC 2 


THE PIRGY OF THERE 39 TRE AOINT S¥ac ThSe ih 
THE REVAIWNER OF TRE FILE CONTAINS YHE ELEVATSONS In reey 


THE TEWWAIN IS THEN MUTPLY TO PAPER TAPE POR INPUT TO PACER ARE. 


ny 


“Te 


; 
ate ae. 


WHITE (60 32) NAF, TITLE, N, ACL2) » CNS, CV CNBOT HS), 10S, 18) NSO ehh) | 


RY’) FORMAT(LML 26%, 242, At, 7, 3HR, BAB, * 
S*ANUMFER OF POINTS wo Massa 
S*GINTERVAL IN INCHES ww %,14,/, 
S"2PGLNT*, 20x," ELEVATIONS Ih FEE *, 
S (79 T5,1PF6,2)) 

8urau, 

SuMeay, 

Cu &S Jai,Nn 
Suwe@iim @ VOT) 
an*geS M2 @ VCCI) awe 
¥(7)82 *& YQ) 

38 Churltuave ‘ 
RMSASHRIC(SUM2 = BUPANe/ND/N) & 320 B 
whT1E(6,36) FMS ’ 

3A Fiber aT(*Arrs w §,610,4," IRCHES%) 
Call FXFOCCS, tater) 

WT 10 1b, 47) PAPAS e A(T2),°0¥(01),T85N) 

GQ Fubra (280, A107, 215,60 0, 3RF 5,99) 
irae Vo ? 96 


ADDED YO THE DISK ON "OM, 


5 ie eg A 


vb. 


a 


Bork PING eL 


i Ager SUBROUTINE PETRICY,' CELYA,OISP) 
| apres ¢ 
: Opad ¢: NETWE: AS Str Tras 
, vars C (1) STARTING AND ENOING POIN(S EQUAL ZEA 
'  apee C (9) MAXIPUM POINT « DSP 
’ aee7 C (3) MINIMUM POINT » DISP z 
‘ Agen Cc 
ARDS Rest ¥ (9) DELTS,OLSP,¢,XK,8 
Poe Te¥eoeR wn, B8Tart, J,0F81,8 
voir C 
seol2 ST4k Ta} 
“VQ1s 1® CONT TAUE 
Ae14 C SCAN UATYL A LIPE IS FOUND THAT UFITS! THE Nata 
aie ie) Bev (STebhy) 
A i 3 Jar, START o2,¢) 
ANT? mBeCY(CITI@VISTARYS) £ (CI eSTART) HDEL TA) 
Agte FaCvtay, / SQRTC1. + Man2) 
AGL Co 2A JeSTARTe1, Tat 
Aaa. xB(TJ=STARTIWOELTA 
uel OEARS(YCII=Ch wXB) w FACT 
Ar2e TF (0 GE. DISP) GO TO 3? 
Aged 2¢ CONTIALE 
Aavza C FCust Cre H 
0225 GO 10 4A 
Vues RY CONTIALE 
a At~ mAVE LONGEST LINE POSSIBLE 
ages ay COMTINE 
ayes PaC¥CTI@«V¥CSTART) 4 (CISSTART) wOELTA) 
wed? FacTes. / SORTCI. + Mea?) 
ogo U2 66 JaeSTaRT,leos 
agso YCJIECVOI) = (Mw (JeSTARYS MOELYA ¢ B29 w FACT 
Acad ou CUP TINUE 
asa START SI 
es) IF (STAR) QLY. Nw1) GO TO &@ , 
AgGK Yinjay, : 
ay 37 IF (START oLTe N) YCNed log, 
We34 RETUQY 
aeady Ern 
veae Ets 


wute LIS) END wwre 
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Pare 


FORTRAN 


‘ 


COMPTL ER SEV. LEV, Jae 


INTEGER FILE C3Y, INT, NUM, REC! (10) 
INTFGER TAPE,PLAT,NAME WRIT, CLOS,OISK 
SCALED FRACTION REC (1H) 

EQUIVALENCE CINT,RECI(C1)9, (NUM, RECLOIM)) 


RATA TAPE/4/, PLAT/22/, NAME/23/, WRIT/16/, CLOS/18/. NISK/I21/ 


Call GMOND (PLAT,OISK,1) 


1A CONTIANTE 
READ (TAPE, 20) FILE, NUM, INT 
OO FORNAT(DAP,A1,/,215) 
TYPR 2H, FILE,NUM,INT 
CALL AMONG (NAME, SISK,FILE,3,1A€) 
Call AMON (WRIT, OTS* ,INT, NUM) 
JEN ehUM / 1A 
Or 49 Je,,JEND 
REANCTAPE, 34) RFC 
AE FORMAT(10S5) . 
CALL IMONDCWRIT,OISK, REC (1), REC(1B)) 
4¢ CONTYAUE Gy 5 
TEND BNR @ JEND MLA . 
YF (TEM FO, AY GO YN 5a 
¢ GET LAST RFCORO 
REAQCTAPE, 30) CRECCI), Te1,IENO) 
CALL SMONOCHRIF,OFSK,RECIS),REC(18)) 
$a CONTINVE 
CALL YMGND(CLOS,OTSK) 
Gr To 1A 
Fun 
PRONGRaM SIZE = "978A 
; CONSTANT T445LE 
257 $5 ant Y 4 
: 260A $17 art 13 
i 241 $eO ncT XY 12 
j 2490 $214 not T a 
: PRAGRAY ALLOCATION 
’ 
\ NAMF aOR NAME ADR NAME ADR NAME ER 
. VT .uyvan 267 v I cies 6 v I OTSK ? AT PILE 1a 
. vt I 2466 VY TENN 265 V YT INT 13 VId 1964 
a VY JEM Pa3 VY NAMR 4 Vv Y NUM 24 V I PLAT f 
AS REC 28 A Ty REC) 13 V Y TAPE 2 V I WRIT i 
; STATEME’'T NUMBER LACATIONS 
; 44. 65 a2 152 ,3e 163, 40 251 ae 
" EXTewh A! DRFERENCTES 98 oa 
: eae 
i ae AS 


PASE 3 Ca-THIS IS THE CONTROL PROGRAM FOR DIGITAL TERRAIN STMULAYS@R am | 


FOE TRAY COMPILER GEV, LEV, Jan 


J\TFEGER KEY, NCHAN,NUMC(A), INIT (9) 
LOGICAL TERM 
RFalL IMY(R) 
COMMON /PARAM/ TERRA, SPACE ,NRUN, DELAY, NPHASE,LAST, INIT 
INTEGER TESRA(H,ZA) ,SPACE(B8, A) ,NRUNCB, 8S DELAY (8,8) pNPHASE (8) © 
LOGICAL LasT (ad 
PATA KEY/OS 
Cer INITIALYZE @#halL OG CONSOLE 
CALL QSHYIM (TERR, 584) 
CALL GSSPCTERR) 
FALL BSCLR (TERR) 
TYPE 7 
7 FORMAT(SEARFNTFER LOGICAL UNIT FOR INPUTS ',/, 
8 BAHKEYROARD = 2',/, 
§ BOHHSPT REANFR o 4!,/) 
ACCEPT 10, KEY 
(A FORMATC(ITS) — 
TEQMakKEY FQ, 2 t 
CALL INPHTCTERM,KEV,INT, NIM) 
TYPE 20 
Om FORMAT(AOHTERRAYN YNPUT COMPLETE!, 7) 
CALL RAYS (TERM, KEY, NCRAN, INT, NUM) 
TYPE 3H 
SA FORMAT(ROHKAAY SETUP COMPLETE !, 7/3 
RITE (HOF, AAS 
40 FPSPMAT(SMHITRAY PHASE TERRAIN SPACE RUNS CHARNEL OELAY') 
ISTARTs} 
wat 
AD CONTYNIHE 
NO GA YweYSTANT,A 
TP (( AST (ID) GO TO 7H 


Ko CONTINUE 
TYFe 6) 

At FHPMAT(AAHERROR 611) 
Caut. EXIT 

7a CONTINUE 


NP aN PHASE (33 
CO 84 Jay ,nP 
HRITE LIAR, BR) KS,TERRAC Ie T) «SPACECI4T) NRUACI,I), 
$ CLADFLAV CI LY ,LeISTART, 2) 
Ag FORMAT(OW, TI AX, TL ,GK, TL, 4X, TS ,1N TK, 4X 51 ,4X,T5, 
% TOA, 38%,54 54,7895 
84 CONTINUE 
JF (CY FG. NCHANS FO TO 9A 
TSTawTael o 1 
Kak «¢ { 
GO Th BA 
Qe CANT TINIE 
Cait GSSOPLJERR) 
CALL ASRUN TERR) 
Celi. -RUNCNCHAN, NIM) 
FAA 99 ve 


PaGE 1 


KORTKRAR FOwPT) &R REV, LFV. Jaa 


SUMOONT TNE INPUT (TERM, KEY, INT, NIM) 
CeaRFANS TS TERSAIN ATA 
LOURTCAL TERY 
INTEGER KEY, NUM(B) 
RFAL INY (8) 
SC4LFD FRACTION SAMPLE (5y9n)° 
TNTFRER AREGIN(B), START 
INTEGER NTFRWA NAME(S), NPLAT,PLAT, POST, READ, DISK, RECL CIA) 
FRuyvslE nce CINTT,RECI (1), CNUML,RECE (IQ) ) 
COMMON /NATA/S BEGIN, SAMPLE 
NATA START s1/ 
NaTa PLAT/9O2/, POSY/2t~/, REAN/I7/, NISKsI2;5/ 
YF (TERM) TYPE Q 
9 FORMAT (BYKENTER NUMBER OF TERRAING ILS) 
REANCKEY,1@) NYEWRA 
1% FARMAT(TS) 
CAL! MANO TELAT, OFSK,1) 
DN 4A Yet,NTERRA ia 
TF (TERM) TYPE 49, 1 


“40 FORMAT(BMNENTER FILE NAME POR TERRAIN NUMBERS!,19,/9 
REANO(KEY, 20) NAME 
2a FORWAT(3a2) 


CALL OMAND(PCSY,DISK,NAME) 
CALL QMOND (REAN, DISK, INTI,NUMI) 
THT CT) SFLOAYCINTL) 12, /828R, 
MM CTY ONUMY 
REGINEI) sSTART 
nA Me JEL NUIMG 10 
Celi. GMOND (READ, DISK, SAMPLE( START) , SAMPLE (START #69) 
STAKTSSYARYT + 190 
3a COM TINUE 
4? CONTINIE 
Be TAK 
ENS 


PROGRAM SI7F #2 1349 


CONSTANY Tate 


M17 $5 ocr yy 3 
320 NOs ocr y 12 
3? #392 ney R Buauna 
322 Orr 4 e 
32% $397 ACT R Slama ‘ ; 
394 Orr 15 : 
PROGRAM ALL OACATIAN 
100 . 

NAME 4anR NAME ANR NAME ADR 
VY ,UtMAe 3497 VT otis 334 Vv 1 OY8k 18 oe a 4 
AO kT 2 Vv I INT! 14 V.eI J 33@ VT KEY 


Pace 4 


FOR T9EE fe. Cyl bw Ry LEV Lae 


TE eae ae AAYSCTFR EY HCHAN,YNT, NUM) 
CeeSFT& He A wAYTAlIM OF FouY Asys 
ChMHCR PPAR AM, TERE A, SPACE, * RUN, PE] AY, NPHASE LAST, ITALY 
Uiarcatl Tees 
TTR GR H KEY, MCHAN, BILE CHY) ,INTT (BR) 
; IS TEGER TAT KAN NAYI FS, TNIMSNTFERRA, RUNS 
RR AL TAY SL Ls Dy ANGLE, MPH MILES, PT,V, WAST 
JS Tesko TERVACAGHS SPACE (RYH) ERUNC RYE) DEL OY (A,B), NPHASE (BD 
MORITA Pasay 
PFAL TNT RY 
PATA PTASLALSO/, YSTAPTSL/ 
TANCTIMSTE CTY 4 Crs) 
CROTCH, AL ESNWT (HWW? © beep) 
ACE EOIN AYEATEND (OPP (AY, A) 
1% Foam eT (ig) 
Jr (TFR™) TyPE aa 
oo FOQMATCRAHENTER “MBER OF RAYS3!,/) 
PLAN TICEY, 194) NRAYS ; 
AO Fae Kw, uBAVS 
TF (TERM) TYPE AG, x 
an FOPMATCBOREL TER KUMAR R CF AXLES FOR BAY', 12,4h8,/9 
RFED (KE Y,I@) NAXLES i. 
MCh ANETSTARY + DwNAXLES = 4 
NN AA \Ve1,NAaAXLFES 
TF (THOM) TYPE an, » 
ae FORMAT CAARENTIFR CISTANCE (FY) FOR AMLEV,YO,1Rigs) 
PEANTKEY, 8a 4 
fal hOWMATCRAS9) 
Limjen 4 SoBe, 
te CONT UML 
i ed On 
r(t)ea, 
Te CTEYH) TYPE Ze, « 
79 FORMATIROHENTER NUMREW CF PHASES FOR BAY',T2,tH3,/) 
FFAQCKEY,1@) “TEYRA 
MO (QRH Jet mTERPA 
TF (TERM) TYPE aa . 
ne FORMAT (ROPENTER TEPPATN, ANGLE (DEGREES) , MPH, MILES, 8',/9 
SFAL (KEY, 9G) TMIIM, ANGLE, MPH, MILES 
On FORMATLT 1, 1X, 3E6,2) 
THETAREANGLE @ PI /18M, 
Nelt « TANCTHFYA) 
VEMPH/3AAD, / {Ona, 
WATTEAINT (TNT (Talay / (VeCNS (THETA) = (8) 
VUESSYETYCMYTLESS/(vesaTT) / FLOATCNUMCTNUM)S » 253 
OC LAM TeISTART,NCHAN 
TERRACJsYISTNUM 
SPACE (JS, II aYFIxX (wath) 
NRUN CS, TI eRUANS 
PHASE (T)EMTEPRA 
LAST(Y) es, FALSF, 
Ima CONTINUE 101 
(ASTINCHAN) BL TRIER, 
TRITCNCREANIETSTART 


es 


et eS aE 


aa ae aa aoe : 


bi 


PAGE 


Te YSTART 
At 111A Neq_sNAKLES 
MELAY OI, TISYFINCL OM) /V ¢ 8) + 2880 
DELAY (Ts YeVOTEIX CCL ONIEOD/V © 68) © 28088 
Jaf o 2 : 
CORTYNUE 
TRE TAMALOSICTNT (TAU), vewanyT) 
ANAL EMOTNHEYTA «9 FRU LSPI 
NaTATCTNUM) / COS(TRETA) 
“WOME £ (KAITSL AND SSAC.) 
MILESBFLOATCOUNS) wFLOATCAUMCTNUM) ) «dD 
TYPE 419, “,JS,ANGLE, MPH, MILFS 
FPORATCRARCOMPIITFO FOR BAY',Y9, BAH, PRASE',I2,/, 
FA,2,9H DEGREES ,F6.2,5K MPH , F672, 6W MILES) 


0 CONTINUE 
JSTARTENCHAN 4 J 
{3 CoH TISukr 


RE THAN 


een 


PRIGRAKR ST7JE & '§526 


CONSTANT TARLE 


1433 
44a? 
yas 
yana 
7435 
14356 
1437 
144” 
44al 
14482 
1443 
§ 444 
1445 
1448 
4447 
4 ann 
AR 
14h? 


gk gery ts 
$1? ncy ¥ 2 
$j4e net ¥ +e 
56% ecy R 512070 
orrT 1$ 
eae ocT Ra 
ony aan 
e613 Oct R §8ane 
Act 10 
eG OCyY &@ 7208 
CCT 44 
$008 OCT P 7H4NA 
mcr 12 
$437 OCT R 4@2ae 
ACT a 
“©7979 OCT | a 
R77 CCT lL. saneae 
$794 OCT T 372a 
APNRMAY ALLOCATION 
NAME ane NAME AOR NAME aor NAKRE 
V I guserse F810 vo Ll gudare 1823 VR ,UgAaR 1584 VR .USBOs 
Vou ,Upoa2 $847 VR ,u200% 3824 VR ,U2@it 16S VR .U2d19 
VR ,UQrF9R 4897 V & PPS 1824 ¥v R ANGLE 1467 V RO 
v7 s425N 6 R INT 3 Y Y ¥START 7 Vz J 
Vv 1K 1658 VI «EY 1 ARL 18 VALI 
Vow FILES 3373 yo R MPR 1a73 Vv IN 1488 V I NAYLRS 
VT RBAYS 184% V J NCRARK ? V I NTERRA 1461 a ¥ NUM 
vok § Vv Y RUNS 1462 102 Vv L TERM B VR THETA 


, OT 


PAF 4 Cue TRUTS NYTGITAL TERRAINS 


" FORTOAN COMPTILEwW SEV, LEV. JAA 


Me 
i 
he 
" 


SUBROUTINE PUM C(NCRAN, NUM) 
TNTEGER NCHAN,NUM(8),INIT(8) 
SCALE FRACTION SAMPLE (C5420) 
JATEGFO RERIN(AY 
d INTEGFR TEPRACB,A), SPACE (CK, 8),NRIUN (8,8), OFILAY (8,8) ,NPHASE (8) 
: LOGICAL L4&ST(B) : 
¥ THYERER NTFRRA, START, NOW, OLE 
: T§TEGER PRASF(R), TIME (A), INDEX (HS, NSEC(B),FIRST(A), STOP (8), 
; % Pri“ (@) 
\ COMMON /PSRPAM/ TERRA, SPACE, NRUN, DELAY, NPHASE,LAST, INIT 
‘4 Coeenn /DATA/S BERLIN, SAMPLE 
: DATA START/S1 BAAS, PHASE/B4IS 
fi &5 Yat,NCHAN 
: \TEOPABSTFRRA(4,T) 
SEC CT) @SPACFE (1,7) 
L@eBEGIN(NTERRA) 
INPEK (VV Sb 
¢ FIRST CT) wl 
: STOPC(T}&. + NUIMCNYERRAY 
ONIM ETP enPUN(y,T9 
TIME (TI RENOFLAY (1,2) 
68 CONTINUE 
CALL GwTIMYCSTARY, JERP) 
CooMaIN UNO 
7A CONYTITAMIE 
CALL ORTIMT (MOW, TERR) 
NLOSSTARYaNOW 
: CALL OwTIMT (START, TERR) 
: CeeShAN CHANNELS 
\ PO RP Jey,nCwan 
t TIME (CYIMTIME(T) © OLN 
‘ IF (YTIME(T) .GT. 7) GO TO 8a 
! CeaQuTPit NISPLACEMENT 
JaINOFY(T) 
CELL AWIPAS(SAMPLE(IJ), Tei, TERR) 
JeJ+ 
: TIME (YJ BNSEC(T) 
TE CF ONE, STCP(I}) GO TO 78 
CewF IN TSHEN nia 
; QNUMCI) SPNUM CT) @] 
IE (PMN CT) .6T, 7) GA TO 74 
Cee TINISHEN PHASE 
JePHASF IT) +4 
16 (3 .GT. NPHASE(Y)) Jat 
NTERDABTERRACY, TY) 
NSEC (1) e8PACE(JS,T1) 
LeGEGIN(NTERRA) 
FIRST(TISL 
STOP (CTI al + MUM(NTERRA) 
ReRUMOC TI aANARIING I, TD) 
PRASF(I)aJ 
Ter C.MOT, LASTCI}) GO Yo 73 
Cae AVY HAS FINISHED PHASE 103 
KYBINITCY) 
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CwemONTPUTS NYSTVAL TERRAINS 


09 79 Kaky,7 
TIME (K) MOELAY(J,K) 


79 CONTINUE 
GO YN %4 
Ce-SFT FOP LONG wart 
7a CONTINUE 
TIME (Po esnane 
74 CONYIUNUE 
JaFTARSTCT) 
74% CONTINIIF 
INDEX (YT) wy 
ay CONTINUE 
Gh TA Ze 
Erp 
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CONSTANT TAALE 


1134 4&8 

{135 $Sa 
T1396 Raha“ 
{137 $545 


PROGRAM ALLOCATION 


NAME 


.Utaan 


Y 
1 
’ 
Y “SFC 
T 
T 


PrPrP<cee< < 


ang 


1] 44 
1147 
1154 
33 
3 
1% 


COMMAN ALL ACATYON 


"tt53 


PPAQAR /ALLOCAT TAM 


VAME 


A Y TERRA 
4 1 NPWASE 


ane 


? 
avo 


NAME AMR 
Vt .Uteia 4148 
aoT INnFy 23 
VEL 1146 
V TNTERRA 144 
A YT Pu 63 


LENGTH w 'a3a 


NAME ADR 


SPACE aa 


A 7 
AL LAST aia 


ADATA /ALLOCATION LEACTH « tyonqa 


MAME 


AT SEIN 


AnR 


i“ 


NAMF Ana 


4 3 SAMPLE ja 


STATEMENT wiiMape LOCATIONS 
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NAME 
4 I FIRST 
‘ I NCHAN 
: T svary 


NAME 


A Y WNRUN 
4 1 INtt 


NAME 


ADR 


43 
115A 


ao-~ 2 


ADR 


280 
4p 


AUR 


NAME 
AY O€Lay 


NAME 


‘aOR 
rT | 


APPENDIX B 
VEHICLE GEOMETRY 
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Ee eRe ER Soe SPRL as MOR ae 
i a re ARNG cee Oe eT OY. Se a ae Tes ae ip 
nc EA PLEA PR Re Ag eS RTT — +: > 


TRUCK CROSSING 


WASHBOARD 
tan @ = — a D*L, tan @ 
cos @ ® =e ee Ws Lhos Q 
Ty we W Ty s I 
Vv a V cos 0 
AND | 
na? ; Ta" 7 | 
L + 
TL2 2 i Tpe s La 9] 
’ v 
Tq = 43 Tp3 * 63_+ 0D 
L3 ° R3 " 2 


U 


( @ - Angle of traversal in degrees (0° + @ < 90°) 
V - Velocity in miles per millisecond 

I - Sample interval in miles 

Ty ~ Time between samples in milliseconds 


Tsy - Initial delays in milliseconds 
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APPENDIX C 
OPERATING PROCEDURES 


fark 20 donee re ELL. Siren ot eed de 25 ae mess Ceo 


OPERATING PROCEDURES 


I. Hewlett-Packard 21168 


A. DETR - Performs detrending and wheel trajectory calculations 
(see section 3 and Appendix A) 


1. The progran first asks for parameters necessary for graphing. 
When the title is asked for, the user should supply blanks 
(the file name is used as the title). 


2. After plotting the terrain profile, the program performs the 


necessary calculations and plots the resulting wheel displace- 
ments. 


3. By entering YES after the wheel trajectory is plotted, a 
listing and paper tape of the wheel profile are produced. 


It. PACER 100 
A. LOAD - Loads paper tapes to disk. 


1. The program continuously loads paper tape without operator 
intervention. 


2. The computer should be halted at the end of the tape. 
B. SAMP - Runs digital simulation (see section 2 and Appendix A) 
1. The EAI 580 should be in digital mode. The program will 
initiate POT-SET and Logic CLEAR. Hydraulics should be 
running but with inputs at zero. 


2. In the terrain input phase enter the one to six character 
name of each file. 


3. The distance for the first axle is its width. The distances 


ae the following axles are the number of feet from the first 
axle. 
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The format for TERRAIN, ANGLE, MPH, MILES is (11, 1X, 366.2). 
Free field input may be used if al] commas and decimal points 
are supplied (see example). 


ERE S 21.) * aiid 
> 
s 


SAEs 
mn 


The computer will go to OPERATE and logic RUN when the bay 
setup is complete. The hydraulics should be brought slowly 
in to verify proper operation. 


ee eee 


a 
e 


To terminate, shut down hydraulics and then HALT computer. 


soe le dl 3 ge 


Sa cite ee 


my 
mu 


t 
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CAAMPLE Ut COMPUTER RUN 


#L,SAMPZ, 21 
Sie? CI Pl 


L 

#G,1000 

ENTER LOGICAL UNIT FOR INPUT: 

KEYBOARD = 2 

Fue READER = 4 

ENTER NUMBER OF TERRAINS: 

ENTER FILE NAME FOR TERRAIN NUMBER: 1 
R005 

RDOS DA Pl 

ENTER FILE NAME FOR TERRAIN NUMBER: 2 
ROD06 

RDO6 DA P] 

ENTER FILE NAME FOR TERRAIN NUMBER: 3 
RDO7 

RD07 DA P] 

TERRAIN INPUT COMPLETE 

ENTER NUMBER OF BAYS: 

Stee NUMBER OF AXLES FOR BAY 1: 

evIER DISTANCE(FT) FOR AXLE 1: 

°* 

ENTER DISTANCE(FY) FOR AXLE 2: 

} ow) 

nies DISTANCE(FT) FOR AXLE 3: 

17.5, 

ENTER NUMBER OF PHASES FOR BAY 1: 

3 

ENTER TERRAIN,ANGLE (DEGREES) ,MPH MILES, : 
3,45., 
COMPUTED” FOR BAY 1, PHASE 1 

44.75 DEGREES 10.00 MPH .25 MILES 
ENTER TERRAIN, ANGLE(DEGREES) ,MPH,MILES,: 
2,45.,20.,.25 

COMPUTED FOR BAY 1, PHASE 2 

44.75 DEGREES 20.00 MPH .23 MILES 
ENTER TERRAIN ,ANGLE(DEGREES) ,MPH MILES, : 
1,45.,30.,.25 

COMPUTED FOR BAY 1, PHASE 


3 
44.75 DEGREES 30.00 MPH 25 MILES 
BAY SETUP COMPLETE 
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Fah oe ss TERR ATA SBP CF bls Ss CHARNEL DFLAY 

1 : a G4 3 4 Qr-an 
9477 
AAR 
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9nape 
227¢ 
2442 
PEED 
PRC? 
phase 


AED HV 
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ase 
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SAMPLE OF COMPUTER PRINTOUT 


| : SPACE: Milliseconds between data. points 
RUNS: Number of times terrain is traversed 
, DELAY: Inftial delay in milliseconds 
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